Noninvasive room temperature instrument to measure liver iron stores 
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1 Introduction 

The work of Brittenham [1], Fischer [2], and others 
has shown that biomagnetic susceptibility meas¬ 
urements are a uniquely effective, noninvasive 
method to monitor liver iron in patients with iron 
overload diseases. However, existing instruments, 
based on superconducting quantum interference 
devices (SQUIDs), are too expensive for wide¬ 
spread clinical use. We are developing a prototype 
biomagnetic susceptometer based on room- 
temperature magnetic sensors. This new instrument 
can potentially make noninvasive liver iron meas¬ 
urements available to patients worldwide. 

In developing this room-temperature biosuscep- 
tometer, we face two main technical challenges. 
First, we need to make sure that the noise in the 
magnetic sensing system is low enough to measure 
the very weak magnetic fields produced by the 
magnetic susceptibility response of liver iron. Sec¬ 
ond, we need to configure the applied-field coils 
and magnetic sensors so as to maximize the signal 
from the liver, while minimizing errors due to the 
response of other body tissues. 

In previous experiments, we achieved an 
instrumental noise level equivalent to a liver-iron 
concentration of 30 |ig Fe /g of wet tissue, in 
magnetic susceptibility measurements on a liver¬ 
sized phantom [3]. This noise level, comparable to 
that of existing SQUID instruments [1,2], is small 
compared with the normal liver iron concentration 
of 200-500 pg/g (wet), and much smaller than other 
sources of error typically encountered in liver iron 
measurements [1,2]. 

However, these initial results were achieved using 
an applied magnetic field that reached well out 
from the sensor unit, so that our magnetic suscepti¬ 
bility measurements sampled a relatively large vol¬ 
ume of material. We were concerned that, in meas¬ 
urements on human subjects, this large sampled 
volume would include too large a contribution from 
the lung, increasing the errors in the liver iron 
measurement. 

We have now redesigned the sensor unit so as to 
minimize the error due the lung, without increasing 
the errors due to instrumental noise. In order to 


optimize the sensor design, we used a numerical 
body model to estimate the magnetic susceptibility 
responses from each of the tissues surrounding the 
liver, as a function of the geometry of the applied- 
field coil and magnetic sensors. In addition we im¬ 
proved the magnetic sensors so that we can make 
magnetic susceptibility measurements with lower 
noise, using a smaller sampled volume. This paper 
presents some of the results from our numerical 
modeling and sensor noise improvements. 

2 Methods 

2.1 Biomagnetic measurement technique 

In biomagnetic susceptometry, we apply a magnetic 
field to the patient's abdomen, and measure the 
change in magnetic field produced by the magnetic 
susceptibility response of the body. Using a water 
bag to fill the space between the sensing instrument 
and the patient's body, we eliminate the large back¬ 
ground signal due to the contrast in susceptibility 
between the body tissues and air. We thus measure 
a magnetic field response proportional to the differ¬ 
ence in susceptibility between the body tissues and 
water. When the sensing instrument is properly 
positioned above the liver, this magnetic suscepti¬ 
bility anomaly is mostly due to the paramagnetic 
response of ferritin in the liver [1,2], 

A major challenge in this measurement is that the 
response of the liver iron is at least 10 7 times 
smaller than the magnetic field that we apply. To 
measure such a tiny response, we need to overcome 
several problems. We need, first, to make the ap¬ 
plied field large enough so that the liver iron signal 
is above sensor noise, without excessive heating in 
the coils that produce the applied field. We also 
need to eliminate the noise caused by fluctuations 
of the current in the applied-field coils. Finally, we 
need to eliminate errors due to temperature drifts in 
the applied-field coils and magnetic sensors. These 
problems are relatively minor in a SQUID instru¬ 
ment, where the noise of the SQUID sensors is ex¬ 
tremely low, large and extremely stable magnetic 
fields can be produced by superconducting coils, 
and the thermal expansion of the coil structure is 
essentially frozen at liquid helium temperatures. 




All of these problems are significant, however, in a 
biomagnetic susceptometer that works at room 
temperature. We overcome these problems by (1) 
canceling the applied field at the magnetic sensor, 
(2) modulating the sample-sensor distance and (3) 
using an oscillatory (ac) magnetic field so that the 
magnetic susceptibility measurement is made in a 
frequency regime where the environmental and 
sensor noise is low. With these techniques, a room- 
temperature susceptometer can achieve noise levels 
comparable to that of existing SQUID systems [3], 

2.2 Optimization of sensor geometry 

Previous work using SQUID sensors shows that the 
errors in actual liver-iron measurements are much 
larger than the noise in the sensing instrument it¬ 
self. For example, in the work of Brittenham et al., 
the standard deviation repeated liver iron measure¬ 
ments on the same subject is in the range of 200 
pg/g of wet tissue [1], It is not known precisely 
what sources contribute to this uncertainty, in what 
proportions. However, it is likely that important 
contributions come from the magnetic susceptibility 
of the lung and the tissues overlying the liver. 

The lung has an important influence because it is a 
large volume of tissue, right next to the liver, with a 
magnetic susceptibility much different from that of 
water. This large susceptibility anomaly can pro¬ 
duce a large magnetic susceptibility signal, that is 
superimposed on the liver iron signal, and varies 
strongly with the position of the sensing unit over 
the patient's abdomen. The overlying tissues are 
important because the magnetic susceptibility 
measurement responds most strongly to material 
close to the magnetic sensors. The response of the 
overlying tissues can be approximately canceled 
out by putting a water bag between the sensing unit 
and the patient's body. However, the response is so 
great that even small variations in the susceptibility 
of the overlying tissues can potentially cause large 
errors in the liver iron measurement. Such suscep¬ 
tibility variations could arise from variations in the 
fat content of the overlying tissue, since fat has a 
slightly less diamagnetic susceptibility than water. 
Variations in the overlying-tissue signal could also 
arise, potentially, from spatial variations in suscep¬ 
tibility due to structures such as the ribs. 

These two sources of error, the lung and the overly¬ 
ing tissue, create a trade-off in the design of the 
applied field coils and magnetic sensors. In order 
to minimize the error due to the lung, it is desirable 
to configure the field coils and sensors so that the 
magnetic susceptibility measurement is sensitive 
only to a restricted body region, which is centered 


over the liver and does not spread out too much, 
laterally, into the lung. However, in restricting the 
lateral extent of the sensed region, we also tend to 
produce a sensitivity profile that dies away quickly 
with depth, thus reducing the liver iron signal rela¬ 
tive to that of the overlying tissue. 

In order to find a best compromise between these 
two effects, we calculated the magnetic susceptibil¬ 
ity response for the liver, lung, abdominal tissue 
and other organs for various configurations of the 
applied field coil and magnetic sensors. These cal¬ 
culations used a digitized model of a human body 
provided to us by Prof. George Zubal of Yale Uni¬ 
versity. This digitized model divides the body into 
243 equally thick slices from head to toe. Each 
slice is further divided into cubical voxels 4 mm on 
a side. Each voxel is assigned an index identifying 
the type of tissue present at that location. Based on 
this tissue classification, we assign each voxel an 
appropriate magnetic susceptibility from Ref. [4], 
To calculate the magnetic susceptibility response of 
a given body region, we first used the law of Biot 
and Savart to calculate the applied field at each 
voxel location, for a given configuration of the ap¬ 
plied-field coils. We then calculated the measured 
magnetic field due to the susceptibility response of 
the tissues, using the formula 

Aj(r ) ( 3B a (r ) • (r - r ) B a (r.) 
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Here, the summation runs through each of the vox¬ 
els in the body model. Av is the volume of each 
voxel, r is the position vector at the center of the 


voxel, r s is the position at which the magnetic field 
is to be measured, B a (r y ) is the applied magnetic 
field and Aj(r.) is the difference in magnetic sus¬ 
ceptibility between the body tissues and water. This 
approximation, which assumes that the total mag¬ 
netic field can be replaced by the applied magnetic 
field at each point in the sample, is valid whenever 
the magnetic susceptibility is small, as it is for the 
diamagnetic materials that make up most body tis¬ 
sue. 

Using this formula, we calculated the magnetic sus¬ 
ceptibility response from each body region, and 
expressed it in terms of the resulting error in the 
measurement of the liver iron concentration. For 
example, to estimate the error due to the lung, we 
first evaluated Eq. (1), including only those voxels 
classified as lung tissue. We assigned each voxel a 
magnetic susceptibility equal to the difference in 



susceptibility between lung tissue (-4. lxlO' 6 SI 
units) and water (-9.05xl0' 6 SI units) [4], We then 
evaluated Eq. (1) again, using only those voxels 
classified as liver, using a susceptibility of 1.6xl0' 6 
SI units, which corresponds to liver iron at a con¬ 
centration of 1 mg/g (wet) [5]. Finally, we divided 
the lung signal by the liver-iron signal to express 
the error due to the lung, in mg of iron per gram of 
wet liver tissue. 

We performed similar calculations for the lung, the 
gas in the intestines and the abdominal tissues over- 
lying the liver. To estimate the signal due to the 
overlying abdominal tissues, we actually calculated 
the magnetic susceptibility response for the entire 
body, relying on the fact that the greatest contribu¬ 
tion to the whole-body signal comes from the tis¬ 
sues immediately below the sensing unit. We did 
these calculations for a number of different designs 
of the applied-field coils and magnetic sensors. 
The following section presents some of the results 
of these calculations. 

3 Results 

Fig. 1 shows the error due to the lung, as a function 
of sensor position, for one of our own sensor de¬ 
signs and for the SQUID sensor now in use in the 
laboratory of Prof. Gary Brittenham at Columbia 
University [1], The error due to the lung is ex¬ 
pressed in mg Fe/mg of wet tissue. The sensor po¬ 
sition coordinate increases as the sensor moves 
down the body (inferiorly), starting from 0 at a 
point 30 mm below the uppermost (superior) tip of 
the liver. In this calculation, the front surface of the 
sensor is tangent to the body surface, at a point ro¬ 
tated 45 degrees toward the right from the midplane 
of the body. 

For the existing SQUID sensor (solid curve), our 
calculations indicate that the lung produces a sub¬ 
stantial error in the liver iron measurement. This 
error is comparable to the normal liver iron concen¬ 
tration, and varies strongly with the position of the 
sensor. In contrast, our own sensor design (broken 
curve) produces a much smaller lung error, which 
varies much less strongly with position. 

First, all of our calculations are based on one body 
model, which incorporates the body dimensions of 
a single individual. Further work is needed to en¬ 
sure that the low lung and abdominal-tissue errors 
in our sensor design are not due to an accidental 
cancellation, dependent on the specific geometry of 
this one person's body. 
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Figure 1: Error due to the lung, in mg of iron per g 
of wet liver tissue, for our own room-temperature 
sensor design and an existing SQUID sensor. 

As shown in Fig. 2, our design also reduces the er¬ 
ror due to the overlying abdominal tissue, relative 
to the existing SQUID design. Generally, we ex¬ 
pect the errors due to the overlying abdominal to 
increase, as we restrict the sampled region in order 
to minimize the error due to the lung. Hence, it is a 
significant gain to reduce both the lung error and 
the overlying-tissue error in the same sensor design. 


Error Due to 1% Variation in Overlying Tissue 
Susceptibility: 

SQUID vs. New Sensor Design 



Figure 2: Errors in liver iron measurement due to 
overlying abdominal tissue, assuming a tissue sus¬ 
ceptibility 1% different from that of water. 

Fig. 3 shows the error in the liver iron measurement 
due to gas in the large intestine. This error is also 
smaller for our sensor design than for the existing 
SQUID system, although the maximum error in 
either case is less than 100 pg Fe/g of wet tissue. 
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Figure 3: Error due to gas in the large intestine, for 
the existing SQUID and our new sensor design. 

In interpreting these results, we need to bear in 
mind the following points: 

Second, to estimate the error due to overlying tis¬ 
sue, we have assumed arbitrarily that the suscepti¬ 
bility of the overlying tissue differs by 1% from 
that of water. At present, the actual variability of 
the tissue susceptibility is not well known. Some 
information is available from the work of Fischer et 
al., who have used a two-parameter fit to estimate 
the susceptibilities for the liver and overlying tissue 
[6], Fischer et al., found overlying-tissue sucepti- 
bilities of-8.6xl0' 6 to -8.8xl0' 6 (SI). Further work 
along these lines may help to characterize more 
completely the variability of the susceptibilty re¬ 
sponse of overlying tissues. This information will 
help us to find the right balance between minimiz¬ 
ing the error due to the lung and minimizing the 
error due to overlying tissues. 

Our calculations, like the two-parameter model of 
Fischer et al., assume that the overlying tissue can 
be characterized by a single value of magnetic sus¬ 
ceptibility. This assumption is certainly question¬ 
able, since the abdominal wall contains several 
types of tissue - muscle, fat and bone — all of 
which have different magnetic susceptibilities. The 
ribs are a particular concern, because they are com¬ 
pact structures, close to the sensing instrument, 
with a susceptibility significantly different from 
that of most other body tissues. The ribs could po¬ 
tentially produce a susceptibility signal that varies 
strongly as the sensor is moved over the surface of 
the body. To evaluate this potential problem, we 
calculated the susceptibility signal as a function of 
sensor position, including only the response of 
bony tissue. We found that the bone signal was 


small compared with that due to liver iron, and var¬ 
ied by only a few percent as the sensor was scanned 
inferiorly down the body. This result suggests that 
the ribs do not produce significant errors in the liver 
iron measurement. 

Finally, in the specific body model used for these 
calculations, the distance from the liver to the outer 
body surface is approximately 2.5 - 3 cm. This dis¬ 
tance is not atypical for the general adult popula¬ 
tion, but is somewhat larger than is typical for the 
thalassemia major patients on whom most biomag- 
netic liver iron measurements are now made. This 
relatively large liver-surface distance significantly 
reduces the signal due to the liver, so that the sig¬ 
nals from the lung and abdominal tissues produce 
proportionately greater errors in the liver-iron 
measurement. In some sense, therefore, our calcu¬ 
lations exaggerate these errors. However, as bio- 
magnetic susceptibility is applied to a wider range 
of patients, including those with hemochromatosis, 
it will be important to design instruments that can 
make accurate measurements in patients with liver- 
surface distances typical of the general population. 
It is in this area that advances in sensor geometry 
will be particularly important. 

As we adjust the sensor geometry in order to mini¬ 
mize errors due to the lung and overlying tissue, we 
vary the depth to which the susceptibility measure¬ 
ment samples the body. This change in sample 
depth affects the magnitude of the liver iron signal, 
and thus changes the effect of sensor noise on the 
liver iron measurement. However, by combining 
model calculations with experimental results, we 
can estimate the effect of sensor noise in a sensor 
system of any desired geometry. Using an im¬ 
proved magnetic sensor, we recently measured the 
susceptibility of a liver-sized phantom, with a root- 
mean-square (rms) instrumental noise equivalent to 
16 fig of iron per g of wet tissue. For the sensor 
geometry in Figs. 2 and 3, we calculate that the 
liver response would be approximately half of that 
in our recent experiments. This reduced liver signal 
would increase the error due to sensor noise by ap¬ 
proximately a factor of two. 

4 Conclusions and future work 

Our calculations and experimental results indicate 
that the room-temperature biosusceptometer will 
have an instrumental noise level of approximately 
30 pg of iron per g of wet tissue. This noise level is 
comparable to that of existing SQUID systems, and 
small compared to other sources of error in the liver 
iron measurement. Furthermore, our model calcu¬ 
lations indicate that this low noise can be achieved 




in a realistic sensor geometry, for which the errors 
due to the lung and overlying tissue are smaller 
than those in the existing SQUID system. 

We are now building a room-temperature biomag- 
netic susceptometer based on the design calcula¬ 
tions presented here. During the next few months, 
we plan to test this instrument on human subjects, 
in the laboratory of Prof. Gary Brittenham at Co¬ 
lumbia University. Our first experiments will in¬ 
vestigate various sources of error, in order to opti¬ 
mize the geometry of the sensor unit. We will then 
compare liver iron values obtained using both the 
room-temperature instrument and the existing 
SQUID system, in a population of patients with 
normal liver iron and iron overload. 
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